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ABSTRACT 

The dipolar magnetic field structure of the neutron star is modified by the plasma originating in the 
pulsar magnetopshere. In the simplest case of a stationary axisymmetric force-free magnetosphere, 
a self-consistent description of the fields and currents is given by the well-known pulsar equation. 
Here we revise the commonly used boundary conditions of the problem in order to incorporate the 
plasma-producing gaps and to provide a framework for a truly self-consistent treatment of the pul- 
sar magnetosphere. The generalized multipolar solution of the pulsar equation is found, which, as 
compared to the customary split monopole solution, is suggested to better represent the character of 
the dipolar force-free field at large distances. In particular, the outer gap location entirely inside the 
light cylinder implies that beyond the light cylinder the null and critical lines should be aligned and 
go parallel to the equator at a certain altitude. Our scheme of the pulsar force-free magnetosphere, 
which will hopefully be followed by extensive analytic and numerical studies, may have numerous 
implications for different fields of pulsar research. 

Subject headings: pulsars: general — stars: magnetic field — stars: neutron — plasmas — magneto- 
hydrodynamics 



1. INTRODUCTION 

Pulsars are the rotating magnetized neutron stars 
(|Hewish et alj [l968) with the rotation periods of 1CU 3 — 
1 s and the surface magnetic field strengths of 10 9 — 
10 12 G. The induction electric field can extract charged 
particles from the stellar sur face and supply the pulsa r 
magnetosphere with plasma ()Goldreich fc Jul ian 1969). 
In the tube of open magnetic field lines, the longitu- 
dinal electric field accelerates the particles to the en- 
ergies of 10 1 2 — 10 13 eV, enabling them to launch 
pair cascades (jRuderman fc Sutherland! [1975). The re- 
sultant secondary electron-positron plasma screens the 
accelerating field all over the tube, except for several 
small regions of pair production called gaps (see Fi g. [T]i, 
Ruderman fc Sutherlandlll975l: IAronslll983l: ICheng et all 



1986). The closed field lines are equipotential, and the 
primary plasma they bear screens the accelerating longi- 
tudinal electric field everywhere. 

The presence of the plasma modifies the original dipo- 
lar structure of the pulsar magnetic field. The prob- 
lem of a self-consistent description of fields and cur- 
rents in the pulsar magnetosphere was formulated in 
the form of a well-known pulsar equation (I Michell 1 1 9 73at 
iScharlemann fc WagoneHll973t lOkamotolll974D . The ba- 
sic model is that of a stationary axisymmetric force- 
free dipole, where the magnetic and rotational axes are 
aligned, the electromagnetic forces are balanced and the 
particle inretia is ignored. An exact analytical solution 
of the pulsar equation was fou nd only for the magnetic 
monopole located at the origin ()Michellli97 3a). in which 
case the field remains monopolar everywhere. In the case 
of a magnetic dipole, the pulsar equation was solved in 
terms of se ries for several specific forms of the current dis- 
tribution (|Michellll973bl IBeskin et all 119831 ). However, 
these solutions are valid only inside the light cylinder 
and cannot be smoothly continued to the infinity. 
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A certain progress was ach ieved in the pioneering work 
of iContopoulos et al.l (|1999f) . which was devoted to the 
numerical treatment of the pulsar equation. As a part 
of the numerical procedure, the current distribution was 
constructed so as to provide a smooth behaviour of the 
magnetic fi eld lines all over t h e spa ce. The numerical 
solution of IContopoulos et al.l (|T999i) was su bsequently 
reproduced by means o f other algorithm s dGruzinov 
1 2005] iKomissarovl 120061 : iMcKinnevI 120061: iSpitkovskv 
120061 ) and generalized in a number of aspects. In 
particular, the non-axisymmetric case was addressed 
( Spit kovskvll200 6: Kala potharakos fc Contopqulosl [20091 : 
iBai fc Spitkovskvll2010tlKalapotharakos et al.ll2012l ). the 
polar gap potential drop and the resultant diffcr- 
ential rotation of the magnetosphere were included 
(|Contopoulosll2005t lTimokhinll2006L 120071 ). the non-ideal 
magnetosphere with a finite c onductivity wa s considered 
(jKalapotharakos et al.l l2012t iLi et all 120121 ) . All these 
studies confirmed the formal validity of the original re- 
sults obtained by IContopoulos et all (I1999D . However, 
it was soon recog nized (jTimokhinl l2006| )~that the polar 
gap can hardly be responsible for the numerically simu- 
lated current distribution, in which the reverse current 
partially flows on the open field lines controlled by the 
gap. This questions the very existence of the stationary 
force-free configuration in pulsars and challenges both 
the polar gap theories and the pulsar magnetosphere 
models. The difficulty of including the polar gap into 
the global structure of the pulsar magnetosphere stimu- 
lated numerous studies of the non-stationary polar gap 
(iLevinson et all 120051: ILuo fc Melrosi 120081: iLvubarskvl 
120091 : lTimokhmll2010| ). In the present paper, we develop 
another approach to solving this difficulty. Namely, we 
revise the physical model underlying the numerical sim- 
ulations of the pulsar force-free magnetosphere. 

We turn to the analytical treatment of the pulsar equa- 
tion and find the multipolar solution, which generalises 
the monopolar one. Based on this result, we suggest a 
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new model of the dipolar force-free magnetosphere of a 
pulsar, which includes all the magnetospheric gaps into 
both the boundary conditions and the current circuit 
configuration. Section 2 contains the basics of the ex- 
isting model. In Sect. 3, we search for the multipolar 
solutions of the pulsar equation. An improved geometri- 
cal model of a force-free dipole beyond the light cylinder 
is constructed in Sect. 4. The force-free magnetosphere, 
which properly includes the magnetosphereic gaps, is de- 
scribed in Sect. 5. Our results are briefly summarized 
and discussed in Sect. 6. 

2. BASIC EQUATIONS 

A starting point for the studies of pulsar electrodynam- 
ics is the model of a stationary axisymmetric force-free 
magnetosphere. It is convenient to choose the cylindrical 
coordinate system (p, <fi, z) with the axis along the pul- 
sar axis. Then the pulsar equation reads (iMichell fl 973a; 
IScharlemann fc Wagonerlll973t [Okamotoll 19741) ~ 
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where the dimensionless functions f{p, z) and A(f) are 
proportional respectively to the magnetic flux and elec- 
tric current through the circle of a radius p centered at 
the magnetic axis at an altitude z above the origin. Both 
functions are unknown, and the problem lies in finding 
the current distribution which makes the self-consistent 
magnetic field obey the physically meaningful boundary 
conditions. 

Strictly speaking, an exact form of the boundary con- 
ditions is not known in advance as well. The set of as- 
sumptions usually taken in the numerical simulations is 
schematically presented in Fig. [T]d. At the neutron star 
surface, the magnetic field is dipolar, / = p 2 /(p 2 + z 2 ) 3 / 2 ; 
at infinity, the field lines become radial, / = f{p/z); 
beyond the light cylinder, the separatrix between the 
open and closed field line regions goes along the equa- 
tor, df/dp = 0; the closed field lines cross the equa- 
tor perpendicularly, df /dz — 0. The current function 
Ad A I df is zero on the closed field lines and at the mag- 
netic axis; the return current flows along the separatrix 
and on the neighbouring open field lines. It should be 
noted that this standard set of assumptions does not al- 
low for the presence of the plasma-producing gaps in the 
pulsar magnetopshere, and the magnetospheric structure 
derived on its basis cannot be completely self-consistent. 
The question of a proper inclusion of the gaps into the 
global magnetospheric structure will be addressed below. 

3. EXACT SOLUTIONS OF THE PULSAR EQUATION 

An exact analytical solution of the pulsar equation ([I} 
is known only for the case of a magnetic monopole lo- 
cated at the center of a neutron star. In order to gener- 
alize this result, we searched for the multipolar solutions, 
which satisfy the relation 
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Fig. 1. — Schematic structure of pulsar magnetosphere: a - the 
case of an axisymmetric vacuum rotating dipole; b - the exist- 
ing force-free model; c - the model suggested in the paper. Solid 
lines represent the magnetic field lines; dashed lines show the light 
cylinder boundary, where the particle rotational velocity equals the 
speed of light; dotted lines delineate the null line, where the charge 
density necessary to screen the accelerating electric field changes 
the sign; thick lines correspond the the critical line, where the cur- 
rent density becomes zero; the direction of current is shown by 
arrows. 

Based on Eq. ([3]) , the magnetic flux function can be pre- 
sented as f(p,z) = f[p^{z)\, where £(z) is an arbitrary 
function. Using this in Eq. © yields 



and 



f = fo 



f = f Q p 2 (z-a), A = 2f 
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The solution 
It represents 
axis by a dis- 



where /o and a are arbitrary constants, 
given by Eq. (|4]) is of particular interest, 
the field of a monopole offset along the z- 
tance a (see Fig. EJd), with the current function being 
identical to that of a monopole centered at the origin. 
From the mathematical point of view, such a generaliza- 
tion of the monopole solution could be directly expected 
from the form of Eq. (fTJ), since it does not depend on the 
z-coordinate explicitly and a shift along this coordinate 
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Fig. 2. — Subsidiary force-free configurations: a - the model of 
a centered split monopole. The level lines of the flux function / 
are shown together with the level values, b - the field of an offset 
monopole according to Eq. Q; /o = 1, a = 1. c — the split-offset 
monopole model; /o = 1, a = 1. The critical field lines and the 
coincident null lines are shown by thick lines. 

changes nothing. The implications of the offset monopole 
solution for the magnetosphere of a force-free dipole are 
discussed in Sect. 0] 
In the spherical coordinate system (r, 9, ip), where r = 

y/ p 2 + z 2 and 9 — atan(p/z), the solution (0} can be 
presented as an infinite sum of the centered magnetic 
multipoles, 
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1 — cos 9 
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yT - 2(r/a) cos9 + r 2 /a 2 
K [P fe _i (cos 9) - cos 9P k (cos 6)] , - > 1, 



where Pfc(cos(9) are the Legendre polynomials. This re- 
sult can be directly obtained from the pulsar equation 
written in the spherical coordinate system by expanding 
the solution into a series in 1 /r k and taking into account 
Eq. @. Note also that although the pulsar equation is 
non-linear, the polynomial function in /, standing on the 
right-hand side, is also a series in 1/r . At r — >• oo, the 
solution coincides with that for the centered monopole. 



But at finite distances the critical line, / = /o, goes par- 
allel to the equator at an altitude a (cf. Fig. [2)3). 

4. GEOMETRICAL MODEL OF A FORCE-FREE DIPOLE 
BEYOND THE LIGHT CYLINDER 

Although the assumption of the neutron star field in 
the form of a monopole is non-realistic, it is gener- 
ally believed that the version of a split monopole (see 
Fig. [2k), where the solution in the upper half-plane is 
symmetrically continued beyond the magnetic equator, 
well mimics the str ucture of the force-free dipole at in- 
finity (|Michellll99lH . However, at finite distances it does 
not allow for the consequences of the essentially dipolar 
features, such as the presence of the magnetospheric gaps 
and the closed field line region. This will be amended in 
our geometrical scheme based on the offset monopole. 

The analogy with an offset monopole enables to ac- 
count for the presence of the outer gap in the pul- 
sar magnetosphere. This gap forms at the intersection 
of th e open field lines with the null line (Chen g~et al.l 
Il986f). along which the mag netospheric charge density 
(e.g. IContopoulos et al.|[l999l) 
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equals zero. It is natural to assume that the outer gap 
is located entirely within the light cylinder. In order to 
provide a proper outer gap location, the null line should 
go along a certain field line all the way beyond the light 
cylinder. Furthermore, if one take that the outer gap 
produces the reverse current necessary to close the pulsar 
circuit, then the null line should coincide with the critical 
line, at which the current function AdA/df equals zero 
(see Fig. [TJ;). Finally, these two lines should be parallel 
to the equator (cf. Eq. ([5])), similarly to the case of an 
offset monopole (see Fig. [2d). Therefore we suggest that 
the force-free field of a dipole at large distances is better 
represented by the split-offset monopole scheme based 
on the two symmetrically offset monopoles of opposite 
polarity (see Fig. [5t). 

Further generalization of the monopolar case can be 
made from analysing the puls ar equation which i ncludes 
differential rotation (see, e.g.. IContopoulosll2005D . 
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where SI = £l(f) is the angular velocity of the magneto- 
sphere rotation allowing for the potential drop across the 
open magnetic field lines. Integration of Eq. ([6]) for the 
case of a monopole yields 



A = 0/ (2 - ///„) . 



(7) 



It is interesting to examine the split-offset monopole con- 
figuration for different pairs (A, SI) obeying Eq. ([7]). As 
the magnetic and electric field strengths read 
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one can see that along the equator 

B 2 -E 2 = - . f ° (9) 
(a 2 + p 2 ) 2 y ' 

and along the horizontal null (critical) lines, f = fo, 

B 2 -E 2 = fjp\ (10) 

i.e. in all these cases B 2 — E 2 = B 2 > (where B p is the 
poloidal held) and the force-free approximation is valid 
for any pair (A, CI). Moreover, as the right-hand sides 
of Eqs. (|9|)- ([10|) arc independent of z, the equilibrium 
condition for the thre e above mentioned horizontal lines , 
d{B 2 - E 2 )/dz = (IOkamotol[T97l ILvubarskvlll990h . 
is also fulfilled for any (A, ft). Thus, the four regions 
bounded by the three horizontal lines in Fig. [2£ may 
have different (A, ft). 

The split-monopole scheme contains from one to three 
current sheets located along the null (critical) lines and 
equator. In the simplest case, A = fl = in the equato- 
rial region between the lines f = fo and A ^ between 
the magnetic axis and these lines. Then the symmet- 
ric current sheets along the null (critical) lines close the 
current circuit in each hemisphere. Given that in the 
equatorial region A ^ as well, the regions on both 
sides of the null (critical) lines may join without a cur- 
rent sheet, in which case the return current flows along 
the equator and the equatorial field lines entering the 
equatorial current sheet. Note that a similar structure of 
the equatorial region beyond the light cylinder is charac- 
te ristic of the dipola r force- free magnetosphere simulated 
in lGruzinovl (|2011al fbT). The only distinction is that in the 
dipolar case B 2 — E 2 = along the equator. It should be 
kept in mind that the realistic magnetosphere of a pul- 
sar may bear not only quantitative but also qualitative 
distinctions from the simplistic monopolc-based scheme. 
In particular, a simpler magnetic field structure in the 
equatorial region, which may result from reconncctions, 
is not excluded. Nevertheless, the split-offset monopole 
scheme suggested above is believed to give insight into 
the general structure of the pulsar magnetosphere allow- 
ing for the coexisting gaps. 

5. NEW MODEL OF DIPOLAR FORCE-FREE 
MAGNETOPSHERE 

Based on the split-offset monopole model, the global 
structure of the stationary axisymmetric dipolar force- 
free magnetosphere of a pulsar seems to look as follows 
(see Fig. [T]:). The critical field line, which becomes par- 
allel to the equator beyond the light cylinder, divides 
the overall open field line region into the two parts. The 
upper part, where the field lines are inclined to the mag- 
netic axis at an acute angle all the way up to infinity, 
is controlled by the polar gap. The lower part, where 
the field lines ultimately enter the equator, is controlled 
by the outer gap. The two gaps are adjusted by the slot 
gap located between them. Note that a similar configura- 
tion of the gaps acting at different bund les of open mag- 
netic fi eld lines was recently obtained in lYuki fc Shibatal 
<|2012f ) by means of numerical simulations of the plasma 
particle motions in the pulsar magnetosphere. 

The current circuit looks as follows. The direct cur- 
rent flows along the field lines controlled by the polar 
and slot gaps and returns to the star in the equatorial 



current sheet, which stretches up to the light cylinder, 
along the equatorial filed lines and through the outer 
gap. Alternatively, beyond the light cylinder the current 
may return in the current sheet along the null (critical) 
line and inside the light cylinder through the outer gap. 
In this case, the longitudinal current flowing through the 
outer and slot gaps changes along the magnetic lines. 
Such a change may be compensated for by the trans- 
field currents in the vicinity of the light cylinder, so that 
continuity equation may still be satisfied. 

Presumably, the outer gap is not infinitesimally nar- 
row, in which case the force-free consideration is no 
longer appropriate. Moreover, the equatorial region be- 
yond the light cylinder is suggestive of the dissipation 
processes, which should be studied out of the framework 
of the force-free approximation in conjunction with the 
outer gap physics. It is important to note that there may 
be a physical relation between the upper and lower hemi- 
spheres through the outer-gap-controlled field lines. We 
also anticipate that the non-stationary dissipation pro- 
cesses in the equatorial region may be responsible for the 
magnetar-like activity of pulsars. 

6. DISCUSSION AND CONCLUSIONS 

We have found a new exact analytic solution of the pul- 
sar equation, which generalizes the well-known monopo- 
lar solution to the case of an offset monopole and presents 
an infinite series over the centered multipoles. We ar- 
gue that, as compared to the classical split monopole, 
the geometrical model based on the two symmetrically 
offset monopoles of opposite polarity better mimics the 
character of the dipolar force-free field beyond the light 
cylinder. In particular, it allows for the presence of the 
outer gap in the pulsar magnetosphere. In our model, 
the polar and outer gaps control different bundles of open 
magnetic field lines, the two gaps being adjusted by the 
slot gap located between them. Then the direct current 
flows through the polar and slot gaps and returns to the 
neutron star through the outer gap. Generally speak- 
ing, in the outer and slot gaps the longitudinal current 
may change along a field line because of the trans-field 
currents in the vicinity of the light cylinder. 

The main distincti on of our force-free mod el from the 
customary one (e.g.. iContopoulos et al.lH999l ) lies in the 
equatorial magnetic field s tructure. In this respect, our 
model is similar to that of iGruzinovl (|2011 a.b). who was 
the first to suggest the alternative boundary condition 
at the equator based on his simulations in the context of 
strong-field electrodynamics. Of course, those works did 
not incorporate the plasma-producing gaps into the mag- 
netospheric structure, but it is the attempt to allow for 
the finite conductivity that have demonstrated the ne- 
cessity to reconsider the equatorial boundary condition. 

Our schematic model of the pulsar force-free magne- 
tosphere including the magnetospheric gaps is believed 
to be a proper framework for detailed analytic and nu- 
merical studies. In particular, we plan to give an exact 
analytic description in the forthcom ing papers. The first 
paper of the series (|Petrovall2012[ ) concentrates on the 
region close to the magnetic axis, and the main result 
is that the force-free field at the top of the polar gap 
differs from that of a vacuum dipole. This is attributed 
to the action of the transverse current which flows at the 
neutron star surface and closes the pulsar current circuit. 
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Our model may also give insight into a detailed pic- 
ture of the particle flows supporting the necessary global 
currents and into the underlying physics of the mag- 
netospheric gaps. Being brought into correspondence 
with the global magnetospheric structure, geometry and 
physics of the gaps should have weighty implications for 
the properties of both the high-energy emission and the 
secondary plasma produced there. It will become possi- 
ble to finally establish which gap is responsible for the 
observed non-thermal high-energy emission and what is 
the underlying physical mechanism. This is especially 
important in the light of the recent progress in pulsar 
studies at GeV energies (e.s.. lAbdo etatl l2010h . The 
refined secondary plasma characteristics would be use- 
ful for elaborating the theory of radio wave propag ation 
in the magnetosphere (jPetrova fc Lyubarskvl 120001 ) and 
may provide a clue for understanding the pulsar radio 
emission mechanism. 

The model of a stationary force-free magnetosphere 
can be regarded as a starting point for studying the in- 
stabilities in the plasma flow and the resultant diversi- 



form fluctuations in pulsar radio emission over a wide 
range of timescales, including those yet to be discov- 
ered in the present and future pulsa r transient searches 
at the world largest radio telescopes ( Konovalenko et all 
l2010t IStappers et al.ll20ll ILiu et al.lboill) . Analysis of 
the pulsar energy losses related not only to the current 
losses on the open field lines, but also to the energy re- 
lease via reconnecti ons (Lyubarsky 1996) and outbursts 
(|Contopoulosl [2005) in the equatorial region, would un- 
cover the relevant scenarios of the neutron star spin 
evolution, cast ne w light on the magnetar-like activity 
(TMereghettil 12008ft and enable one to construct a physi- 
cally grounded classification of the neutron star observa- 
tional manifestations beyond the fr amework of the clas- 
sical rotating radio pulsar concept (jKaspi 2010). 



I am grateful to Prof. S. Shibata for useful discus- 
sion. The work is partially supported by the grant of the 
President of Ukraine (the project of the State Fund for 
Fundamental Research No. F35/554-2011). 
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